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. Abstract. We present 2D results of simulations of the magnetorotational 

core collapsed supernova. For the first time we obtain strong explosion for the 
core collapsed supernova. In 2D approximation we show that amplification of 
the toroidal magnetic field due to the differential rotation leads to the formation 
of MHD Shockwave, which produces supernova explosion. The amounts of the 
ejected mass 0.1M® and energy ~ 0.5 4- 0.6 • 10 51 ergs can explain the energy 
output for supernova type II or type Ib/c explosions. The shape of the explosion 
is qualitatively depends on the initial configuration of the magnetic field, and 
\ may form strong ejection neat the equatorial plane, or produce mildly collimatcd 

• jets. Our simulation show that during the evolution of the magnetic field the 

magnetorotational instability appears and leads to exponential growth of the 
magnetic field strength. 
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1. Introduction 



CO 

The reliable explanation of the mechanism of core collapsed supernova is still 
open question for the modern astrophysics. We describe the results of 2D nu- 
merical simulation of magnetorotational mechanism for core collapsed super- 
nova. The main idea of the magnetorotational mechan ism for the core collapsed 
supernova was suggested by Bis novatvi-Koganl (|l97Cf l. After the core collapse 



central parts of the star rotate differentially. In the presence of initial poloidal 
magnetic field toroidal component of the magnetic field appears and amplifies 
with time. When the pressure of the magnetic field becomes comparable with 
the gas pressure the compression wave arise and moves along steeply decreasing 
density profile. In a short time this compression wave transforms to the MHD 
fast shock wave and produces supernova expl osion. First 2D simul a tions of the 
magnetorotational mechanism were made by Le Bl anck k, Wilson! (^70) with 
unrealistically large initial magn etic field. 2D s imula tions of the s imila r prob- 
lem were given in the papers of lOhnishil (I1983T ) and ISvmbali^tvl (Il984l) . The 
magnetorotational processes with very large initial magnetic fiel ds (typical for 
magne tars) were recently simulated numerically in the papers of iKotake et alJ 
(2004]) and lYamada fc Sawail (l2004fi . The ID s imulations of the magnetorota - 



tional supernova were presented in the papers bv lBisnovatvi-Kogan et all (|l976t ) 
I Ar deli an etafl (|1979T ). In the last paper the simulations were performed for c 



wide range of the initial values of magnetic field. 2D simulations of the col- 
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Figure 1. Velocity field for t = 1.06s after beginning the evolution of the 
toroidal magnetic field component for the dipole-likc initial magnetic field (left 
plot), velocity field for t = 0.2s after beginning the evolution of the toroidal 
magnetic field component for the quadrupole-\\ke initial magnetic field (right 
plot). 
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Figure 2. Time evolution for of the ejected mass (in relation to Mq) and 
energy for the quadrupole-\ike initial magnetic field, a — 10~ 6 . 



la pse of the magnetized rotating protostellar cloud are presented in the paper 
bv lArdelian et al.l <|200Cf t. 



2. Formulation of the problem 

For the numerical simulation of the magnetorotational supernova we use set 
of MHD equations with selfgravitation and with infinite conductivity in cylin- 
drical coordinates. The e quations of state and neutrino losses are described 
in details in the paper of I Ar deli an et aD ( 20041 ) . Axial symmetry (^p = 0, 

r > 0) and symmetry to the equatorial plane (z > 0) are assumed. For the 
simulations we used specially developed implicit numerical method based on the 
La grangian operator-difference s chem e on triangular grid of variable structure, 
see I Ar deli an Sz Kosmachevskii I (|l995l ). 
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Figure 3. Time evolution of the ejected mass (in relation to Mq) and energy 
for the dipole-\\ke initial magnetic field, a = 1CP 6 . 



3. Dipole-like and quadrupole-like initial magnetic fields 

For the simulations we used initial dipole-like and quadrupole-like initial mag- 
netic fields. Magnetorotational explosion with the initial dipole-like magnetic 
field has the shape of mildly collimated jet (Figure Q left plot). This jet can be 
collimated later when the shock wave pass through the envelope of the star. Ini- 
tial quadrupole-like poloidal magnetic field leads to the explosion which develops 
preferably near equatorial plane (Figure Q right plot). 

In the both cases (dipole and quadrupole initial poloidal magnetic field) 
the amplification of the toroidal magnetic field leads to the formation of the 
compression wave, which moves along steeply decreasing density profile, and 
soon it transforms to the fast MHD shock wave. 

For the quadrupole initial magnetic field the evolution of the magnetoro- 
tational explosion takes place much faster then for the dipole field (Figur^2J 

n. 

The simulations of the magnetorotational explosion for the initial quadru- 
pole and dipole magnetic field lead to the ejection of about 0.6-10 51 ergs of energy 
and about O.14M of mass. The amount of the ejected mass will increase with 
the evolution of the MHD shock. 



4. Magnetorotational instability in the magnetorotational mecha- 
nism 

Let a be the ratio of the initial poloidal magnetic energy of the star and its 
gravitational energy at the moment of the "turning on" magnetic field, a = 
^grav • ^ e nave made simulations for the initial quadrupole-like magnetic field 
for the wide range of the initial values of the magnetic energy (a = 10 -2 -^10~ 12 ). 
The values of the explosion energy and ejected mass are approximately the 
same for all c alculated varian t s. In ID simulations of the magnetorotational 
mechanism by I Ar deli an et al.l (1979) it was shown that evolution time of the 
magnetorotational explosion depends on a in the following way t exp i os ion ~ 
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In 2D simulations the situation as qualitatively differen t. Due to the a p pearing 
and e v olution of the magnetorotational ins tability, see iDungevI (|l958fi , iTavlerl 
(ll97.t ). ISpniitl (|2002l l lAkivama et alJ (I200.TT ). the time from the beginning of the 
amplification of the toroidal component of the magnetic field to the moment of 
explosion is much shorter then in the results of ID simulations. At the Figure|IJ 
the dependence of the time of the explosion (which is calculated as a time from 
the beginning of the evolution of the magnetic field up to the moment of the 
beginning of supernova explosion) from a value is represented. For a >~ 10~ 4 
the 1/y/a law holds approximately, but for smaller values of a the explosion 

time texplosion 

is proportional t exp i osion ~ 
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